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Abstract. The physics of the n phase shift in ferromagnetic Josephson junctions may enable a range of 
applications for spin-electronic devices and quantum computing. We investigate transitions from "0" to 
'V states in Nb/Fe/Nb Josephson junctions by varying the Fe barrier thickness from 0.5 nm to 5.5 nm. 
From magnetic measurements we estimate for Fe a magnetic dead layer of about 1.1 nm. By fitting the 
characteristic voltage oscillations with existing theoretical models we extrapolate an exchange energy of 
256 meV, a Fermi velocity of 1.98 x 10 5 m/s and an electron mean free path of 6.2 nm, in agreement with 
other reported values. From the temperature dependence of the IcRn product we show that its decay rate 
exhibits a nonmonotonic oscillatory behavior with the Fe barrier thickness. 

PACS. 74.50.+r Tunneling phenomena; point contacts, weak links, Josephson effects - 74.25.Sv Critical 
currents - 74.78.Db Low-Tc films - 74. 25. Ha Magnetic properties 

1 Introduction enabled by the development of nano-scale technology [2]. 

In particular, ferromagnetic 7r-j unctions have been pro- 

In recent years a great interest has grown in the imple- P osed as coherent two-state quantum systems 3 j. In a 

mentation of quantum computing [T] , where information is Superconductor/Ferromagnet (S/F) system the supercon- 

stored in two- level quantum systems (qubits). Many the- ducting order parameter penetrates into the F layer. As 

oretical and experimental advances have been achieved a consequence of exchange splitting of the spin-up and 

towards the physical realisation of qubits and, in this re- spin-down electrons in the ferromagnetic sub-bands the 

spect, solid-state implementations appear as good candi- superconducting order parameter oscillates as a function 

dates thanks to their potential scalability which is partly of depth in the F-layer which causes its sign to change pe- 
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riodically |4|5|6j . A positive order parameter corresponds 
to a "0" -state, while a negative order parameter to a "tt"- 
state. To realise a practical quantum device it is neces- 
sary to control these two states; this is easily achieved in 
S/F/S systems incorporating weak [7] and strong [819110] 
F barriers and hence such systems have been identified as 
suitable quantum-electronic devices based on 7r-shift tech- 
nology [11] . These systems exhibit oscillations of the crit- 
ical current with F layer thickness equivalent to a change 
of phase from to 7r-states. The principle goal in this re- 
search field is to understand the physics of superconductor 
based quantum technology and to evaluate potential ma- 
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Fig. 1. (Colour online) Saturation magnetisation vs. Fe barrier 
thickness for three deposition runs. The measurements were 
taken at room temperature. From the linear fit the magnetic 
dead layer is extrapolated. Inset: illustration of our multilayers. 



terial systems, while aiming at reducing the size of the n n . . r nr , ln , ^ ,, -, r 

J o o clean limit condition is fulfilled, hrom the study of the 



employed heterostructures towards the nano-scale. 



temperature dependence of the IqRn product we have 



T ,i • , c , , j c ii iu- i shown that its decay rate exhibits a nonmonotonic oscil- 

In this paper, we present the first study of the thick- J 

i i r , , ... i , .,, ,. . . latory behavior with the Fe barrier thickness. Our analysis 

ness dependence of the critical current oscillations m mo- J J 

bium/iron/niobium (Nb/Fe/Nb) Josephson junctions. Mag- reveals that Fe ' for its sma11 ma S netic dead ^ and hi S h 

, r • r i , , . i exchange energy, is well-suited for the implementation into 

netic measurements as a function of the he layer thickness ° aj 1 

i , \ , , i , ,. , £ ,t. . nano-scale superconductor based quantum electronics. 

itpe) nave been made so that an estimate of the magnetic ^ n 

dead layer of our Fe barrier layer can be obtained. The cur- 
rent vs voltage characteristics have been measured for dif- r» .■ .■ «- ml. ir /ml. 

6 2 Magnetic properties or Nb/Fe/Nb 

ferent thicknesses of the Fe and show periodic transitions 

from 0-states to 7r-states from nodes in the I c RN{tFe) re- Our heterostructures consist of Nb (250 nm) / Fe (0.5 nm 
lation, where I c is the critical current of the device and Rn - 5.5 nm) / Nb (250 nm). To assist subsequent process- 
is the normal-state resistance. From a comparison between ing in a focused ion beam (FIB) microscope, 20 nm of Cu 
experimental data and two theoretical models we have es- was deposited inside the outer Nb electrodes, but 50 nm 
timated the exchange energy and the Fermi velocity of the away from the Fe barrier. We remark that the 20 nm of 
Fe barrier, which are found to be in good agreement with Cu is a thickness smaller than its coherence length, so it is 
the expected values from literature. With a simple model completely proximitised into the Nb, and does not affect 
we estimate the mean free path of Fe verifying that the the transport properties of the Josephson junction. The 
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heterostructures are deposited by d.c. magnetron sputter- 
ing in an Ar plasma at 1.5 Pa on 10 x 5 mm silicon (100) 
substrates coated with a 250 nm thick oxide layer on the 
surface. Our deposition system is equipped with a rotat- 
ing holder which can move under three magnetrons and 
allows the loading of more than twenty substrates in one 
run. In this way our system permits, knowing the relation 
between the deposition rate and the speed of the rotat- 
ing holder, to deposit in-situ and in a single run, different 
samples with various Fe barrier thicknesses. From x-ray 
reflectivity measurements, we have verified the control of 
Fe barrier thickness variation t p e to be within an accuracy 
of ±0.2 nm [12]. 



200 



100- 



-100 



-200 



■L =4.1 nm without microwaves 

Fe 

•t =4.1 nm f= 13.86GHz 

Fe 




-3 



Fig. 2. (Colour online) / vs. V/((f>of), where 4>o = h/2e, for 
a Josephson junction with an Fe barrier thickness of 4.1 nm 
(black dash line) and the same device with 13.86 GHz excita- 
tion of microwaves (red line) at T= 4.2 K. Inset: an illustra- 
tion of a focused ion beam processed device where x ~ 4/xm, 
100 < y < 500 nm and 100 < z < 500 nm. The current path is 
also shown. 



To investigate the magnetic properties of our devices 
we have measured, using a vibrating sample magnetome- 
ter at room temperature, the magnetic moment per unit netic dead layer so that better control and tunability of 
surface area of the films as a function of Fe barrier thick- the Nb phase-state is possible. The issues associated with 
ness. The saturation magnetisation was measured for three a magnetic dead layer are a particular problem with the 
different deposition runs (see Fig. [T]), where for each depo- strong ferromagnetic barriers because of the very thin lay- 
sition we have obtained similar saturation magnetisations ers that have to be deposited in order to exploit or even 
which also confirms our control of the Fe barrier thick- see 7r-shift physics and hence, improved accessibility and 
ness. By extrapolating the least-squares fit of this data knowledge of this dead layer is important if strong fer- 
we have estimated that the Fe has a magnetic dead layer romagnetic materials are to be incorporated into 7r-shift 
to — 1-1 nm. The thickness of the magnetic dead layer is based technologies, 
one of the central points into the realisation of quantum 
devices because it marks the lower limit of the ferromag- 
netic barrier thickness for which it exhibits magnetic mo- 
ment. We remark that the scale of the estimated Fe barrier 

magnetic dead layer is low; however, there is a need for The Josephson junctions were fabricated in three steps: 



3 Realization of nanostructured Josephson 
junctions and electrical characterization 



deposition optimisation to minimise the scale of the mag- (i) patterning of films using optical lithography to de- 
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fine micron-scale tracks and contact pads. Our mask per- 4 Critical current oscillations 

mits 14 devices per sample which allows us to measure 
numerous devices and to derive good estimates of impor- 
tant parameters, like I c Rn] (h) broad beam Ar ion milling 
(3 mAcm~ 2 , 500 V beam) to remove unwanted material 
from around the mask pattern which leaves 4/zm tracks 
for subsequent FIB work; (iii) FIB etching of micron-scale 
tracts[T3] to achieve vertical transport with a device area 
(yz) in the range 0.2 /an 2 to 1.1 /im 2 . We show a labeled 
illustration of a device in the inset of Fig. 



For each sample we have measured, from the I — V char- 
acteristics, at least three junctions from which the critical 
current (I c ) and the normal resistance (Rn) of a device 
is extracted. The I c Rn product is then plotted as a func- 
tion of Fe barrier thickness and reveals I c RN{tFe) decay- 
ing with multiple oscillations in the clean limit up to 5 nm 
(see FigEJ). 

The experimental data have been modeled with the 
following formula [TS]: 



I c Rn = ^c-Riv(io) 



Transport measurements were made in a custom made 
liquid He dip probe with a heating stage and microwave 
antenna fitted. Measurements were made with a standard 
lock-in technique. The critical current I c was found from 
the differential resistance as the point where the differ- 
ential resistance increases above the value for zero bias 
current. Rn was measured using a quasi-dc bias current 
of 3-5 mA, this enabled the nonlinear portion of the I-V 
curves to be neglected, but was not large enough to drive 
the Nb electrodes into a normal state. 

In Fig. [2] we show I c Rn vs. V/<fiof for a device with 
an Fe barrier thickness of ~ 0.8 nm with and without mi- 
crowaves applied, where </>o is the flux quantum given by 
h/2e and / is the applied microwave frequency (/ = 14.01 
GHz). Constant voltage Shapiro steps appear due to the 
synchronisation of the Josephson oscillations on the ap- 
plied excitation [T3] at voltages equal to integer multiples 
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«2 



where t\ is the Fe thickness corresponding to the first min- 
imum in I c Rn, IcRn (to) is the first experimental value 
of I c i?Ar(i_F e ), and £i and £2 are the two fitting param- 
eters. We notice that the experimental data are in good 
agreement with this theoretical model, and from the fit 
we obtain ^ = 3.8 nm and £2 — 0-25 nm and a period 
of the oscillations is t = 1.6 nm. Eq. {1} is a general 
formula that applies to the clean and dirty limits; how- 
ever, in the clean limit £2 = v F h/2E ex . By assuming the 
Fermi velocity of Fe is vf = 1-98 x 10 5 m/s, as reported 
in literature [16], we calculate the exchange energy of the 
Iron: E ex — Kvp/2t ~ 256 meV. 

To confirm that our oscillations are in the clean limit 
(meaning that the considered Fe thickness is always smaller 
than the Fe mean free path), we have modeled our data 
with a second simplified formula which holds only in this 
limit [IT]: 



I sm{2E ex tF /TiVf) I 

of hf/2e = 0.029 meV. IcRn 00 2E ex t F /hvf ' (2) 
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Fig. 4. Dependence of £2/^1 with the inverse magnetic length; 
Fig. 3. (Colour online) I C R N oscillations as a function of the inset . vs L /^ H to estimate the mean free path, L. 



thickness of the iron barrier (dot) with the theoretical fit (line). 

where, in this case, E ex and vf are the two fitting param- 
eters. From the theoretical fit we extrapolate E ex — 256 
meV and vp = 1-98 x 10 5 m/s. 

This model is fitted using exactly the same values for 
the parameters as those obtained from the previous model 
given by Eq. [TJ Both models are thus consistent with one 
another and show an excellent fit to our data, which ex- 
hibit multiple oscillations of I c Rn in a small range of Fe 
barrier thickness. In Fig. Owe show the experimental data 
with the two theoretical fits. The magnetic dead layer 
has not been subtracted, as for T= 4K it is expected to 
be significantly lower than 1.1 nm. This is an important 
feature towards the realisation of integrated nanostruc- 
tured quantum electronic devices based on 7r-phase shift 
in Nb/Fe/Nb Josephson junctions. 

With a simplified model that is obtained solving the 
linearised Eilenberger equations [18j we can estimate the 
mean free path of Fe. The general formula is: 



tanh ■ 



L 



eff 



(3) 



where £ e // is the effective decay length given by = 
^i 1 + ; Co is the Ginzburg-Landau coherence length 
and £h is a complex coherence length. In the clean limit 
1 + LCo 1 » (max{ln(l + L^ 1 ), ln(L^ H 1 )})/2. The so- 
lution of this equation gives: £f 1 = S.Q 1 + L^ 1 , £0 = 
{v F h)/{2irT c k B ) and 

In our case the ratio £2/^1 — 0.06, and so assuming 
L/£o=0.2 we can extrapolate from the graphical solution 
a value of L/£h — 29 (FigHJ). Considering the curve L/& 
vs L/£h (inset in Fig[4| we obtain a value of L/& to be 
~ 25 and hence, we estimate the mean free path of Fe 
to be ~ 6.2 nm. With this analysis we remark that the 
condition that I c RN(dFe) is within the clean limit only is 
fulfilled. 

We have followed the temperature dependence of the 
IcRn product for different thicknesses of the Fe barrier, 
as shown in FigJ51 We note that, for each Fe barrier thick- 
ness, the Josephson junction resistance remains approxi- 
mately constant. On the other hand, the critical current, 
and hence the IcRn product, quickly decreases with in- 
creasing the temperature. It is interesting to notice how 
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barrier thickness on I c Rn- From magnetic measurements 
the Fe barriers were found to have a magnetic dead layer 
of ~ 1.1 nm. From I ~ V curves we have measured the 
Josephson critical current (/ c ) and the normal-state resis- 
tance Rn so that oscillations in I c Rn as a function of the 
Fe barrier thickness could be studied. In agreement with 
two theoretical models we have fitted our data and esti- 
mated the exchange energy of Fe to be 256 mcV and the 
Fermi velocity of Fe to be 1.98 x 10 5 m/s. Furthermore, we 
have applied a simple theoretical model that involves solv- 
ing the linearised Eilenberger equations so that an esti- 
mate of the Fe electron mean free path could be made (6.2 
nm). A value of 6.2 nm confirms that I c RM{tFe) oscilla- 
tions are fully in the clean limit giving rise to a maximised 
I c Rn product in the 7r-state. For different Fe barrier thick- 
nesses we have shown that the IcRn product decreases 
with increasing temperature, and in particular the decay 
rate presents the same oscillatory behavior as the critical 
(4) current, this original feature seems common to strong fer- 
romagnetic 7r junctions [12j . although no explanation for 
this can yet be given. With this work we have shown that 
it is possible to fabricate nanostructured Nb/Fe/Nb tt- 

j unctions with a small magnetic dead layer of 1.1 nm and 
observe oscillations of a T in phase with the oscillations of ^ & Ugh ^ q{ oycr thfi Fe b&rrier thickness 

the IcRn, in other words the smaller IcRn decays to zero mi ; j i r-n • i ^ 

° ' variation. 1 he estimated exchange energy of le is close to 



Fig. 5. (Colour online) Temperature dependence of the IcRn 
for different thicknesses of the Fe layer (logarithmic scale). In- 
set: oscillations of the decay rate cct(To ~ 4.2K) of IcRn, 
Eq. Q, with the temperature [stars; the blue line is a guide to 
the eye], as compared with the oscillations of the IcRn prod- 
uct itself [spheres; the red line is the fit given by Eq. JT}] in a 
window of Fe thicknesses. 

the rate at which the critical current drops to zero exhibits 
a nonmonotonic dependence on the Fe layer thickness. We 
can define the relative decay rate olt of the IcRn product 
with the temperature as in Ref. [8], namely 

dhx[I c R N (T)} 



a T (T ) 



dT 
dI c R N (T) 



T=T n 



dT 



1 



t=t q IcRn(T ) 



calculated at To = 4.2 K. In the inset of Fig[5]we plot ax 
as a function of tp e for five consecutive thicknesses; we 



with increasing T more slowly than the greater IcRn- 



5 Conclusions 



bulk Fc implying that the Fe is clean and S/F roughness 
is minimal. Some interfacial diffusion of Fe into Nb could 
account for the slight suppression of E ex (Fe) and the mag- 



In summary, we have fabricated Nb/Fe/Nb nano-scale Joseph- 



nitude of the magnetic dead layer. We conclude that Fe 



son junctions and investigated the effect of varying the Fe 



barrier S/F/S Josephson junctions are viable structures in 
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the development of superconductor-based quantum elec- 12. Robinson J. W. A., Piano S., Burnell G., Bell C, 

tronic devices. The electrical and magnetic properties of and Blamire M. G., submitted. 

Fe are well understood and are routinely used in the mag- 13- Bell C, Burnell G., Kang D.-J., Hadfield R. H., 

netics industry, therefore Nb/Fe/Nb multilayers can read- Kappers M. J. and Blamire M. G., Nanotechnology 14 



(2003) 630. 

14. Barone A. and Paterno G., Physics and Applications of 
the Josephson effect (John Wiley & Sons, New York) 1982. 

15. Born F., Siegel M., Hollmann E. K., Braak H., Gol- 
ubov A. A., Gusakova D. Yu. and Kupriyanov M. Yu., 

1. Nielsen M. A. and Chuang I. L., Quantum Computa- Phys. Rev. B, 74 (2006) 140501(R). 



ily be used in controllable two-level quantum information 
systems. 
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